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Abstract 

Recently, a set of thermodynamic Bethe ansatz equations is proposed by Dorey, 
Pocklington and Tateo for unitary minimal models perturbed by 0i 2 or (p2,i oper- 
ator. We examine their results in view of the lattice analogues, dilute Al models at 
regime 1 and 2. Taking + (f>i t 2 and M^ t 4 + 02,1 as the simplest examples, we 
will explicitly show that the conjectured TBA equations can be recovered from the 
lattice model in a scaling limit. 
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1 Introduction 



Since the breakthrough in the integrable perturbation theory of CFT Ej, there has 
been a lot of progress in the understanding of 4>i,3 perturbation theory Hj. On the 
other hand, although the remarkable example, the Ising model in a magnetic field , was 
treated in [I], the progress on the 4>x,2 an d 02, 1 perturbed theories has been steady but 
slow. 

The systematic studies on the bootstrap procedure on S matrix have been initiated 
in |Hj and [Hj. The latter approach, based on the scaling q— state Potts field theory, has 
been further elaborated by Dorey et al [7j. Thanks to the Coleman-Thun mechanism, 
they argue that the contributions from spurious poles cancel and conclude the closed set 
of S— matrices for a wide range of parameters. 

The check of the results against a finite size system, however suffers from the non- 
diagonal nature of the scattering process. Due to the lack of a relevant string hypothesis, 
the diagonalization of the transfer matrix is far from trivial. In [Hj, conjectured is a 
set of thermodynamic Bethe ansatz equations (TBA) from the consideration on special 
cases of which they found similarity to the TBA for the sine-Gordon model . Roughly 
speaking, they proposed the TBA by gluing the "breather-kink" part and the "magnon" 
part in which the latter originates from the sine-Gordon model at specific coupling^ ITUj. 
Although the derivation is intuitive, the resultant equations pass many non trivial checks. 

In this report, we shall examine the problem in view of a solvable lattice model . As a 
lattice analogue to M i i+1 + 1)2 , M i+1 £+2 + 02,i we consider the L— state RSOS model 
proposed in JTT| [T2*]. which will be referred to as the dilute model. There are several 
evidences for this correspondence, the central charge jllj. the scaling dimensions of the 
leading perturbation [IH [TH], universal ratios [HJITHIHE! an d so on - 

The question whether it shares the identical TBA to describe its finite temperature 
(size) property has not yet been fully answered. The purpose of this report is to present 
positive evidences for this inquiry. 

There are already few examples to the demonstration of the equivalence. The common 
TBA of the dilute A 3 model at regime 2 and the M 3i4 + 4>i,2 case is firstly proved in [Hj . 
There the most dominant solutions to the Bethe ansatz equation are explicitly identified 
in the form of the "string solution" , which leads to the famous E 8 TBA. In the case 
L = 4,6, corresponding to the E 7 , E 6 case, such explicit identification of string hypothesis 
seems not yet to be completed. 

An alternative approach, based on the quantum transfer matrix (QTM) fEJEI|, has 
been successfully applied to L = 3,4,6 20, 21 . The functional relations among properly 
chosen QTMs play the fundamental role there and it enables to derive TBA without 
knowing the explicit locations of dominant solutions to Bethe ansatz equation . 

For L = 3,4,6 cases, the underlying affine Lie algebraic structure (E 8 , E-j, E 6 , respec- 
tively) provides several clues in the investigation of the functional relations among QTMs. 
The remaining case, which seems to lose a direct connection to affine Lie algebra in general 
(see, however exceptions [S]). It might be thus challenging to the clarify the functional 
relation, and thereby see if the Y system in is actually recovered. In this report, the last 
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"exceptional" case (in terminology of 8.J) M 5 6 for the 0i j2 perturbation, and the first 
exceptional case M 3 4 for the 2 ,i perturbation are focused. 

This paper is organized as follows. In the next section, we give a brief review on the 
dilute Al models and the QTM method. Section 3 is devoted to the discussion on the dilute 
A 5 model at regime 2 which is expected to be a lattice analogue of the M 5j6 + 0i j2 theory. 
Fusion QTMs parameterized by skew Young diagrams are introduced and found to satisfy 
a set of closed functional relations. It will be shown that the conjectured TBA is naturally 
derived in a scaling limit. In case of the dilute Al model, L even, a fundamental role 
seems to be played by a "kink" transfer matrix. As the simplest and the most well-known 
example, we treat + 02,i, corresponding to the Ising model off critical temperature, 
in section 4. We conclude the paper with brief summary and discussion in section 5. 

2 The dilute Al model and the quantum transfer ma- 
trix 

The dilute Al model is proposed in ^T] as an elliptic extension of the Izergin-Korepin 
model [221- The model is of the restricted SOS type with local variables G {1, 2, • • • , L}. 
The variables {a, b} on neighboring sites should satisfy the adjacency condition, [at — bj < 1, 
which is often described by a graph in figO In [H], the RSOS weights, satisfying the Yang- 




Figure 1: An incidence diagram for the dilute A 5 model. The local states corresponding 
to connected nodes can be located to nearest neighbor sites on a square lattice. 

Baxter relation, have been found to be parameterized by the spectral parameter u and the 
elliptic nome q. The crossing parameter A needs to be a function of L for the restriction. 
The model exhibits four different physical regimes depending on parameters, 

• regime 1. < u < 3, A = 

• regime 2. < u < 3, A = 

• regime 3. 3 — ? < u < 

• regime 4. 3 — ^ < w < 

We are interested in regimes 1 and 2. 

The central charge and scaling dimension associated to leading perturbation evaluated 
in [HI HHj suggests, 

• The dilute Al-i model at regime 1 is an underlying lattice theory for Ml,l+i + 02, i 



4(L+1) ' n — ^ 

*(L+2) j ^ o 
4(L+1) ' U — ° 

A - 4(L+1) ' L - 6 
^ = 4(1+1) ' L - 2 " 
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The dilute Al model at regime 2 is an underlying lattice theory for Mj^+i + 



'1,2 



There are also further evidences supporting this correspondence, as mentioned in the 
introduction. 

One can introduce an associated ID quantum system to the above 2D classical model. 
The Hamiltonian 7~Cxd for the former is defined from the row to row transfer matrix 
Trtr(w) of the latter, by 

d 

Kid = £-Q^\ogT RTR (u)\ u=0 . 

We omit the explicit operator form of Ji\D- The parameter e = —1, (1) labels regimes 
1 and 2 (3 and 4). 

The thermodynamics of the ID quantum system is the central issue in the following. 
We apply the method of QTM [THJ HE] to this problem. Leaving details to references, we 
list the only relevant results for the following discussion. 

A fundamental QTM is defined in a staggered manner 



N/2 



b 2 
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In the above, squares represent Boltzmann weights; four indices represent local variables 
and the spectral parameters are specified inside of them. The fictitious dimension N( 
even) , sometimes referred to as the Trotter number, is introduced. It has nothing to do 
with the real system size of the original ID system. The real system size will not appear 
in our discussion as the quantities after taking the thermodynamic limit is of our interest. 

It is vital that two (spectral) parameters u, x exist and that only the latter concerns 
the commutative property of QTMs, [Tqtm( m > x )i ^QTM^j x ')\ = 0- The remaining pa- 
rameter u plays the role in intertwining the finite Trotter number (N) system and the 
finite temperature system (/?) by u = u* = — e-j|. More concretely, the free energy per site 
is represented only by the largest eigenvalue T\(u, x) of Tqtm at x = and u = u*, 

Pf = - lim \ogT 1 {u*,x = 0) 

N— >oo 

The eigenvalue T\(u,x) takes the form 

Ti(u,x) = W(j)(x + -1)4>(X + -%)— I-T+0 + [X--1 )— y— -y- 

2 2 Q(x - p) 2 2 Q(x - p) Q(x + p) 

3. w l.x g(g+jj) 

-l)(f){X 1) ; r— . 

2 >Vy 2 } Q(x + U) 



+ w l (j){x - -l)(j){x - -t)^ (1) 



JV 

Q(x) := ] J h[x — Xj] 

■h[x + (| - u)i]h[x - (| - u)i]\N/2 
h[2i]h[3i] 



h[x] := i9i(iXx), 
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where w = exp(ij^j) (£ = 1 for the largest eigenvalue sector). 

The parameters, {xj} are solutions to "Bethe ansatz equation" (BAE), 

w ^i±j) = Q( x j ~ i)Q(xj + gj) j = i N ( 2 ) 
(f)(xj — i) Q(xj + i)Q(xj — 2i) ' 

From now on we suppress the dependency on u which must be set as u = u* . 

It has been shown in many examples [23] > that the functional relations among "gener- 
alized" (fusion) QTMs offer a way to evaluate the free energy without precise knowledge 
on the locations {xj}. We adopt the same strategy here and shall discuss the functional 
relations realized among fusion QTMs of the dilute Al model below. 

3 QTM associated to skew Young diagrams and 
quantum Jacobi -Trudi formula 

We introduce fusion QTMs associated to Young diagrams. The idea to connect Young 
diagrams and (eigenvalues of) QTM, originated in j2U 123 I2H| is very simple. Let three 
boxes with letters 1,2 and 3 represent the three terms in eigenvalue of the quantum transfer 
matrix ((H), 

T 1 (x)=a+a+a. 

Obviously, the eigenvalue of a fusion QTM can be represented by a summation of products 
of "boxes" with different letters and spectral parameters, over a certain set. The point is 
that the set can be identified with semi-standard Young tableaux (SST) for s/ 3 . We state 
the above situation more precisely. Let /i and A be a pair of Young tableaux satisfying 
fii > Aj, VI We subtract a diagram A from fi, which is called a skew Young diagram /i — A. 
The usual Young diagram is the special case that A is empty, and we will omit A in the 
case hereafter. On each diagrams, the spectral parameter changes +2i from the left box 
to the right and — 2i from the top box to the bottom. We fix the spectral parameter 
associated to the right-top box to be x + + /xi — 2) (or equivalently the spectral 
parameter associated to the left-bottom box to be x — + \i\ — 2) ). Insert a letter 
£ij to the -th box such that the semi-standard condition is satisfied. We denote its 
spectral parameter by Xij. Then the product 



n 

is associated to the Young table. The summation over the tableaux satisfying the semi- 
standard condition then defines 



w= e na... o) 



which is expected to be the eigenvalue of a fusion QTM. 
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The simplest subset of the above is the QTM based on Young diagrams of the rect- 
angular shape. It was shown (201 that for any such member reduces to QTM of 1 x m 
Young diagram, which is related to m— fold symmetric fusion. For later convenience, we 
introduce a renormalized 1 x m fusion QTMs T m (x) by 



J K ' il<l2<-<ir, 



ii 


12 




1m 



The renormalization factor f m , common to tableaux of width m, is given by 

3=1 

Hereafter, for any function f(x), we denote by f(x ± ia) the product f(x + ia)f(x — ia) . 

Then the resultant T m 's are all degree 2N w.r.t. h[x + shift], and have a periodicity 
due to Boltzmann weights; T m (x + Pi) = T m (x), where 

. L _ y , for regime 2 




for regime 1. 
Remarkably, T m (x) enjoys a "duality" 

T f x \ _ \ T 2L-i-rn(x), in = 0, • • • , 2L , for L even 

1 T2L-i-m(x + f i), m = 0, • • • , 2L , for L odd. 

This is deduced from the a 2 nature of the model and special choice of A. We have at 
least checked the validity numerically and assume their validity in this report. The above 
two properties, the periodicity and the duality (JSJ) play the fundamental role in the proof 
of the closed functional relations. 

The real usefulness of T m (x) lies in the fact that any QTM associated to a skew Young 
diagram can be represented in terms of their products. 

Theorem 1. Let T^/xix) be a renormalized T^nix) in (0) by a common factor, 
Xj (x + i(X - Hi + fij + Xj - 2j + 1)). Then the following equality holds. 

%/ x {x) = deti< 3ifc < /1 ' 1 (T w _ Afc _ :;+fc (a; + i(/i' x - \i\ + ^ + X k - j - k + 1))) (6) 

where T m<0 := 0. 

We regard this as a quantum analogue of the Jacobi-Trudi formula. 

By this, apparently T^/x{x) is an analytic function of x due to BAE, and contains the 
quantity of our interest, Ti( special case. The former assertion is not obvious from 

the original definition by the tableaux, but it is trivial from the quantum Jacobi-Trudi 
formula. 

In the same spirit, we introduce A m /a(x), which is analytic under BAE, 

A m/ a(s) := T»/\(x)/ '.{T m >2L{x) -> 0}. 
The pole-free property of A M / A (x) is apparent from (JBJ). 
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4 dilute A5 model at regime 2 as a lattice analogue 

to M 5;6 + <j> 1)2 

For M5.6 + </>i,2, Dorey et al argued the existence of two kinds of particles, 2 kinks and 4 
breathers. For diagonalization of scattering theory, they introduced 2 magnons (massless 
particles), in addition. Explicitly, the Y — system reads, 

3 3 
Y Bl (x±— i) = E B3 (x), Y B[i (x ± —i) = E Bl (x)E B5 (x) 

3 2 1 

Y B5 (x±— i) = E B3 (x)E K2 (x± — ^Zk^Z^x ± — i)E 2 (x) 

3 2 1 

Y ih(x±—i) = E Kl (x± — i)Ei(x± — i)E K2 (x) 

Y K2 {x±^i) = E B5 (x)£W{x), Y Kl (x±U) = E B2 (x)rt l \x) 

Yt(x ± — i) = C 2 (x)Ck 2 (x)C Ki (x), Y 2 (x ± — i) = d(x) 
with 

C a {x) := - | 1 E a (x) := 1 + Y a (x) 

1 + Yjx) 

where a takes one of B\, B3, ■ ■ ■ ,1,2. (Y±, Y 2 are written as Y^\ Y^ in |Hj.) 

We are not starting from Y but rather from the QTM. Corresponding to breathers, 
we introduce "breather" QTM by 



T B3 (x) 



T Bl (x) := T x {x) 

13 12 

= A (8)1) + — i)/(f>{x- —i) 

3 

T Bs (x) := A ( i5 i8i8)/(7J) (x)/0(x± -i) 

11 12 9 

T B7 (x) := A (15i15i8i8)/(14i7j) (x+ — - yi)0(x± -i)) 

T B2 (x) := T 7 (x) 

OK 

T®(x) := A (8i7)/(6) (x + -i) 

then the following relations, referred to as the "breather" T system, hold. 

3 11 12 

T Bl (x±—i) = T (x±— i) + (j)(x - —i)T B . A (x) 

3 8 
t b 3 (x ± — i) = T (x)T (x ± — i) + T Bl (x)T B5 (x) 

3 3 5 

T B 5 {x±—i) = T (x±— i)T (x±— i) + T B3 (x)T Br (x) 

T B2 (x±^) = T (x±^)+T^(x) 



where T (x) = f2{%)- They are originated from the "hidden su(2)" discussed in |2"T] . 

In contrast to the dilute A 3 model (equivalently the E 8 case), the "hidden su(2)" 
structure is not enough to obtain a closed set of functional relations. We then introduce 
another set of functional relations, related to magnons. 

To each nodes on the D4 Dynkin diagram (see fig EI), we associate tm\x),(a = 
1,2,3, 4, m G Z> ) and set (x) = 1. Then we impose a L> 4 related T system among 




Figure 2: The nodes in the D± Dynkin diagram are indexed in the above manner, 
them, in terminology of |2*H] . 

<& ) (a;±^)=*S.i(^S-i(^) + ^ ) (^)n*S ) (^) (7) 

where g[ (x ± jz) — g^\x). In the above by b ~ a, we mean that a and 6 are connected 
on the Dynkin diagram. 

Moreover we set an inhomogeneous truncation, = = and put g^ = g^ = 1. 

Unless one introduces some further condition, the set of functional relations (0) are 
not closed, so can not be solved. Then we demand 



= 4>), t^(x)=t^(x)T B3 (x), 
g?(x) = T (x), g[ 2 \x)=T (x±j). 

The second relation glues the breather T system to the D4 related T system system. 

The above requirements seem to be rather artificial, but they lead to remarkable 
consequences. First, solutions to can be given in terms of QTM appearing in the 
dilute A 5 model as follows. 
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t^(x)=T^\x) tf\x)=T B7 (x) 

tf\x) = A(l2,8,7)/(5,4)(^) 4 4) (x) = A( 5 ,1)(X + ^l) / 4>(x - —I + ^z) 

4 1} (x) = T B5 (x)T B2 (x ± 1) 4 2) (x) = T Bs (x ± ^)T B2 (x) 

4^ ( X ) = T B 5 [X ± Y^) T B 2 (X ± ^) 4^ fa) = T Bs (^) T B 2 0*0 • 

The proof of the above statement is too lengthy to reproduce here. We hope to present 
them with the general discussion of L general 

Second, the following combination of T and t solves the Y — system for M 5 6 + 0i >2 . 

cj>(x - fi)T Bs (x) T Bl (x)T Bs (x) 
XbA x ) - TfTJZ i 11 Y B 3 {X) 



T (x ± n i) T (x)T (x ± n z 



2 



>V = y / x *2 (g) 

**TiW (3) (1) J X 2 W (3) (2) 

t\'{x)g\'{x) t\'{x)g\'{x) 

Yi(x) = ,,, ■ ■ ^r- — nrr- _ ^2(2^ — 



4 1} (a;)4 2) (^)4 4) (^) 4 3) (^)' 

Third, the functions T , t, Y possess "nice" analytic properties. Before stating the 
properties, we need preparations. Note that the Y — system is invariant, for even N, if Y 
is replaced by Y, defined by 

f , r^r^ n s for w < 

\y Bi {x)k{x ± z(l - for m > 

and all other cases, F a = V a . The parameter u' stands for yM. This is due to the definition 
of K, 

K x = Z T7^ a J 

which satisfies k(x ± zt|) = 1. The elliptic nome q' = exp(— r'),r' = 4r is introduced so 
as to respect the periodicity of the Y function on the real direction of x. We denote a 
typical Y equation as 



Y a (x ± ia) = Y\ -b(x ± 27;,) Y\ £c(x ± 



«7c 



Our numerical data indicate that the rhs is analytic and nonzero in the strip Qx G 
[—a, a}. Each element in the lhs also satisfies the same in appropriate strips , i.e., Eb(x) is 
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analytic and nonzero in the strip Qx G [— 7b, 7b], and so on. These remarkable properties 
enable us to solve the coupled algebraic equation, like (jHJ), in the Fourier space ( to be 
precise, its logarithmic derivatives). Then the inverse Fourier transformation leads to the 
coupled integral equations which yield the explicit evaluation of log Y" a (x). 

To make a direct contact with the TBA result, three further steps are needed. First 
take the Trotter limit N — > oo, uN — /3, (e = —1). Second rewrite \ogE b (x) by log£b(a;). 
Third, take a scaling limit. The stepl is executable analytically, which manifests one of the 
advantage of the present approach. The resultant equations no longer have dependency 
on a fictitious N but only depends on the temperature variable, (3. After the step 2, we 
obtain the equations, in the Fourier space, 





(lot 






(1) 




(~ Bl \ 


M 

















V 


'■■ J 




W 




\ '■■ J 



where C Bl — log(l + y^-) and similarly for others. The quantities with hat indicate that 

they are Fourier transformations. M and Kq are asymmetric matrices of which explicit 
forms are omitted here but can be easily obtained from the y system. The only first entry 
has a nonvanishing inhomogeneous term in the rhs. This reflects the fact that only Y Bl 
needs some trivial renormalization so as to have nice analytic properties. By multiplying 
M _1 from the left, the kernel matrix of TBA, M~ l K$ turns out to be symmetric, remark- 
ably. This property is crucial in applying the dilogarithm technique to evaluate the central 
charge. The inhomogeneous term vector 47r / 3M~ 1 ■ *(1, 0, • ■ ■ ) possesses six non vanishing 
elements. 



4vr/3(2 cosh f\k + 1)(2 cosh ±k - 1) 

W) 
8tc[3 cosh j^k 
(2cosh^fc- l)D(k) 

8tt(3 cosh ^k 
(2cosh^fc- l)D(fc)' 

where we denote by ds 1 for the drive term associated to logYgj and so on. A common 
denominator D{k) denotes 

12 10 6 4 

D(k) = 2 cosh — k + 2 cosh — k — 2 cosh — k — 2 cosh — k + 1. 
v ; 14 14 14 14 

We finally perform the step 3. In view of QFT, the bulk quantity is not of direct 
interest, rather the fluctuation is. We introduce yB x {x) = Yb 1 (x + r"), for example, to 
evaluate quantities near the "fermi surface" with r" = y 1 . Then take a limit q — > such 

that rrikR = 2^,^-1 ^ ■ By r we mean the residue of i/D(k) at h — n/3i. Two quantities 



d Bl 



d B , 



d>Ki 



8tt(3 cosh jjk 
(2coshf\k-l)D(k) 
16tt[3 cosh j^k cosh jjk 



id 



(2cosh^fc- l)D(k) 

47T/3 

(2 cosh ^A;- l)D(k) 



i Ba 



l B 2 



d 



K 2 
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M~ l and K seem to carry the information of S matrices; the elements of M 1 K agree 
with the expression described in |.8j in terms of S matrices, under identification x = W/ix 
in the limit q —>■ 0. The matrix M _1 also encodes the information of the mass spectra. 
When taking the inverse Fourier transformation, the nearest zero to the real axis, k — 
of D(k), is relevant in the "scaling" limit as r" tends to be infinity. Applying the Poisson's 
summation formula, we found a most dominant term, 



d Kl (x) 



■t>7i(3r 



2 cos § - 1 



- — 7" 1 

-e 7 cosh —x 



tukR cosh 9, 



for example, where \x = 9. Note that the relation oc q~? is consistent with the scaling 
dimension A 1,2 = g. One similarly verifies that all other drive terms also take the form 
mR cosh 9 and their mass ratio agree with those in [Hj. 
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m Bl = 2m K cos —n 



TUB, 



Am K COS — 7T COS — -7T 



42 
4 

m^„ = 2m^ cos — n 
2 42 



42 



11 3 

fflfl, = 4mft- COS — 7T COS — 7T 

1 

rriB 2 = Iviik COS — 7T 



Thus the TBA of M 5 6 + 12 theory is recovered from the scaling limit of the dilute 
A 5 model at regime 2. 

Once Y is fixed by TBA, we can also evaluate the free energy from 

T x {x ± = T Bl (x ± li) = T (x ± + F Bl (x)). 

It is readily shown that a "fluctuation" part of the free energy / is proportional to 
^2 J2k I "^fc-Rcosh 9 log(l + l/yk)d9, which is the desired expression. 

5 dilute A<i model at regime 1 as a lattice analogue 

to M 3)4 + 2) 1 

We treat another example corresponding to 02, 1 perturbation theory, the simplest and 
most well studied case, the Ising model off critical temperature, M 3 4 + <p 2) i- The model is 
described by a free fermion, thus is rather trivial in a sense. In view of functional relations, 
however, it is not trivial to derive the simplest Y system Y(x ± i|) = 1 (in the present 
normalization of x), from T\{x) in (pQ) which consists of 3 terms. This model is actually 
one of the first examples, which require a more fundamental object than Ti(x), a box, 
which seems to correspond to a fundamental breather B\. 
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We define 

, v Q{x + 2i) Q{x)Q{x + 3i) 

w v 1 Q(x + i) K 1 Q(x + i)Q(x - 1) w 

- M >t^f (10) 

which has a property common to 7i(:c) namely, it is pole-free due to the Bethe ansatz 
equation. 

More importantly, we have functional relations, 

r K (x±~i) = T 1 {x)+T 2 (x) = 2T 1 {x) (11) 

T K {x±h) = T 3 (x)+T (x)-<P(x±^i)(w 3 + ^)=2(<P(x±li)+<P(x±^))(l2) 
Z z w 6 2 z 

The first equalities are directly verified by comparing both sides in forms of the ratio 
of Q functions. The seconds are consequences of the duality. One then reaches a desired 
relation (jTJ| after proper renormalizations. The first equation, (fTTjl seems to suggest tk{x) 
is related to the kink in the theory; the bound state of kink produces a breather. 

In general L = even case, we find that Tk(x) plays the most fundamental role, which 
will be a topic of a separate publication. 

It is a nice exercise to recover from (fTT|) and ()12j1 . the free fermion free energy, in the 
scaling limit. We shall remark the analytic property, supported by numerics, that tk(x) 
being Analytic and Nonzero in the strip Imx G [—§,§], for that purpose. 



6 Summary and discussion 

In this report, we demonstrate explicitly that TBA for M 5j6 + 4>i,2 and M 3i4 + 2 ,i, con- 
jectured by Dorey et al, are realized in the scaling limit of lattice models. The crucial 
idea is to introduce fusion transfer matrices associated to skew Young tableaux and to 
investigate the functional relations among them. The proofs of functional relations are 
rather combinatorial and lengthy, thus omitted due to the lack of space. They will be 
supplemented in the subsequent paper which discusses the TBA behind the dilute Al 
models, L general J2H]- 

There are still many open problems. The explicit identification of string solutions 
would be definitely one of the most important. The complete study on this will shed some 
lights on the way how to proceed for TBA in the case of perturbed non-unitary minimal 
models. We mention the first step in this direction in [30] . 
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